Liquid Drops sliding down an inclined plane 

Inwon Kim*and Antoine Mellet ^ 

February 26, 2013 
(N 

o 

^^ Abstract 

d We investigate a one-dimensional model describing the motion of liquid drops sliding 

^TJ down an inclined plane (the so-called quasi-static approximation model). We prove existence 

^^ and uniqueness of a solution and investigate its long time behavior for both homogeneous 

^f^ and inhomogeneous medium (i.e. constant and non-constant contact angle). We also obtain 

T~H some homogenization results. 

^ 1 Introduction 

^ 1.1 Equilibrium drops and quasi-static approximation 



cd 



Consider a liquid drop lying on an inclined plane. We introduce a coordinate system such that 

^ this plane is the {x, y)-plane, and we assume that the free surface of the drop (the liquid/vapor 

interface) can be described as the graph of a function {x,y) h- > u{x,y). For small droplets, the 

,__l energy of an equilibrium drop can be approximated by 

> „(^,^) 

5^ Jf{u)^a j l\Vu\^ + ^X{u>o}dxdy+ f f Tdzdxdy (1) 

^^ where a denotes the surface tension coefficient and (3{x,y) is the relative adhesion coefficient 

CO between the fluid and the sohd. If the support plane is inclined at an angle a to the horizontal 

^—^ in the x-direction, the gravitational potential T can be written as 

r = /9g(z cosa — xsina) (2) 



X 
J3 



and the Eulcr-Lagrange equation for the minimization of ^ with volume constraint is the fol- 
lowing equation (known as the Young-Laplace equation): 

{—Am = X — UK cos a + XK sin a in |u > 0| 

(3) 
i|Vu|2 = /3 ond{u>0}. 

with K = pg/(J and where A is a Lagrange multiplier associated with the volume constraint. Note 
that the Dirichlet integral J | Vup dx dy is an approximation of the surface tension energy, which 
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classically involves the perimeter / -y/l + |Vup dx dy (see [3j). This perimeter functional leads 
to a mean curvature operator which would be much more delicate to deal with than the Laplace 
operator appearing in (Is]). 

It is well-known (see 3 ) that when (3 is constant and k > 0, then (Is]) has no solution (this 
can be seen by multiplying the first equation in ([3| by u^ and integrating over {u > 0} - it is 
also obvious that any translation of u down the x-axis will decrease the energy ^). 

This means that on a perfectly homogeneous surface, any drop should slide down the inclined 
plane, no matter how small the inclination of the plane a. This is not the case when /3 is a 
function of x. In that case, one expects some drops to stick to the inclined plane, at least for 
small inclination (or small volume). This sticking phenomenon was rigorously investigated by 
Caffarelli-Mellet 2] when /? is assumed to be periodic with period £ <C 1. 

In any case, large drops will always eventually slide down the inclined plane and it is the 
purpose of this paper to investigate the motion of such drops, when /3 is constant (homogeneous 
surface) and when it is not (heterogeneous surface) . There is a number of papers (see for instance 
[U EJ [71 [12] and references therein) in the fluid dynamics literature discussing the behavior of 
sliding drops. Interesting phenomena are observed. In particular, the shape of the drop can 
change drastically when the velocity is increased, and a singularity (corner) develops at the rear 
of the drop when the velocity exceed a certain critical value (and a cusp may form at even higher 
velocity). 

In this paper, we aim at studying such phenomena in the framework of a (relatively) simple 
model for the motion of liquid drops usually referred to as the quasi-static approximation. We 
will see that this model does allow for singularity formation of the type discussed above. 

1.2 The quasi-static approximation model. 

In the quasi-static approximation regime, the speed of the contact line is much slower than the 
capillary relaxation time, so that at each instant there is a balance between gravitational force 
and surface tension (see [BIZ]). The free surface of the drop is thus at equilibrium at all time, 
and can be described by a function u{x,y,t) solution of: 

— Au = X{t) — uKCOsa + XKsina in {u(i) > 0}, (4) 

where the constant X(t) is determined by the volume constraint 

u{x,t) dx — Vq. 

Along the free boundary (or contact line) d{u > 0}, the equilibrium contact angle condition 
is assumed to be satisfied only at some microscopic scale, and the deviation of the "apparent" 
contact angle 6 = | Vm| from the equilibrium value 6e = y/2P is responsible for the motion of the 
contact line. The relation between the velocity of the contact line and the contact angle 6 is not 
a settled issue, and various velocity laws have been proposed (see [S], [Z], [E]) and studied (see 
[4]). Following Blake and Ruschak [T, we assume that the speed is proportional to O"^ — &l. More 
precisely, we assume that the normal velocity V of the contact line d{u > 0} is given by 

y=i|VMp-/3 ona{u>0}. (5) 

This particular velocity law has the advantage of leading to a gradient flow formulation for Q- 
((sl). This gradient flow formulation will be discussed in the next section, but let us already 



stress out that (Is]) leads to an obvious problem when the drop is sliding down: In the rear of 
the drop (where V < 0), the speed |1^| is bounded by f3 while in the front of the drop (where 
V > 0), the speed can increase without bounds (and is expected to increase if the volume of 
the drop, or the inclination increase). This would lead to ever-increasing support of the drop, 
which is incompatible with m (indeed, we will see Q does not admits positive solutions on large 
domains). This problem can be fixed by noticing that ^ should not be expected to hold in the 
degenerate case where Vw — along the free boundary. In the next section, we will see that the 
gradient flow formulation naturally leads to an obstacle problem type formulation of the velocity 
law of the form (see Q) 

-F+^|Vwp~/3, |Vu|l =0 ond{u>0}. (6) 



The free boundary problem Q-Q has been studied by Grunewald and Kim [6 in the par- 
ticular case a = 0. In that case, we always have |Vu| > along the free boundary, and so ^ 
holds (instead of (|6])). In the present paper, we are interested in the case a > 0, for which the 
degeneracy of |Vu| cannot be ruled out. In that case, the existence of solutions (which is the 
main focus in [5]) is not the only interesting issue. Indeed, one would also like to know whether 
a given drop will stick or will slide down the plane. And when it does slide down the plane, one 
would like to characterize the motion of the drop. We will see that there exist some traveling 
wave type solutions which describe the asymptotic speed and profile of any sliding drops. Finally, 
one would like to determine the effects of heterogeneities of the inclined plane on the motion of 
the sliding drops. We refer to HO] for results concerning the motion of liquid droplets on a 
horizontal heterogenous plane (see also [llj for similar results in the framework of Hele-Shaw 
flow). 

Such issues are rather delicate to investigate in full generality. In this paper, we develop the 
analysis of (|4|-(l6| (and answer all the above questions) in the one-dimensional case. In particular, 
we prove the existence and uniqueness of solutions (for general functions /3) and we study their 
long time behavior when /3 is constant, and when ^ is periodic. This one-dimensional model 
could be interpreted as modeling the motion of a fixed volume of fluid initially placed uniformly 
across the inclined plane (ridge of fluid). However, such a configuration is not usually stable (the 
leading edge may become unstable in the span- wise direction as it flows down the inclined plane) . 
We thus prefer to think of this study as a starting point for a more general study involving the 
higher dimensional case. 

As we will see, even this (simpler) case already leads to interesting behaviors. In particular, 
the profile of the drop may degenerate in the rear of the drop (touching the ground tangcntially) . 
Of course, in dimension two and higher, the dynamics of the drops is much more complex since 
topology changes seem unavoidable (splitting and merging of droplets) . 

Finally, let us mention that one of main difficulty in studying this problem (regardless of 
the dimension) is the lack of comparison principle due to the volume constraint (no viscosity 
solutions), since A depends on i in Q. And there is no way to get rid of this volume constraint, 
especially when addressing the long-time behavior of the drops. Indeed if we were to fix A(i) = Xq 
and solve Q-Q, it is easy to check that in most cases the drop would simply vanish in finite 
time. 



1.3 Gradient flow structure 

As mentioned above, the system of equation Q-Q has a natural gradient flow structure: For a 
given function v, we denote by 

^{v) = / -\\7v\'^ + —Kcosa — V XKsinadx dy 
Jr2 2 2 

the energy of the corresponding drop, and for a given Caccioppoli set W cM.'^, wc define 

uw ^ a.rgmm{^{v) ; veH^{W), v>Q, jv^Vo}. (7) 

When a = 0, it is easy to show that the solution of Q with support W and volume Vq is 
nonnegative in W. When a > 0, this may not be the case for large set W. This is the reason 
why we need the additional obstacle condition w > in ([7]). As a consequence, we may not have 
supp uw = W, which is the main difference with the case a = 0. 
We now define the energy of a Caccioppoli set W as: 



£•{]¥) = ,^{uw)+ / (Sdxdy 



w 



Formally, Q-lls]) is the gradient flow for the energy (o over the set of Caccioppoli sets with respect 
to the usual Riemannian structure. Indeed, let us consider a solution W{t) of the gradient flow 

dtW = -diff f^{W) 

and denote by V the normal velocity of W. We can show that V satisfies ([5| as long as W = 
suppuvy: Let f; be a normal velocity field applied to dW and let dii be the variation of uw 
induced by v. Assuming that |-Duvi/| j^ on dW, we get 



di f f S" (W) (v) — / Duw ■ DSu + {u\YKCOsa — xKsma)Sudx + 

Jw Jd 

= / [ — Auw + uwi^cosa ^ xKsma\Sudx + 

Jw Jd 



P+^\Duw\^ 



V dS 



'W 

X Su dx + 



w 



aw 



\Du 



w\ 



p+-\Duw 



dW 
2 



-\Duw\5u + ( /3 + -\Duw\^ ) V 



dS 



vdS 



1, 



l3--\Duwr]idS, 
aw \ ^ 



where in the first equality we used the fact that uw = on dW, and the third inequality is a 
consequence of the following expansion: 



This implies 



Su ~ v\Duw\ + o{v) on dW. 



V^-\Duw?-p. 



(8) 



However, as mentioned above, when a > 0, the solution of y\ may have a support strictly smaller 
than W when W is large (because gravity is pushing u " downward" , in the direction of increasing 



x). In that case, the obstacle free boundary condition yields |VMvy| — along d{uw > 0} \ dW 
and we do not recover ([S]). If we denote by V the normal velocity of supp uw, it is clear that 
whenever supp uw ^ W, we have V < V. Since this can only happen when |Vun/| = 0, we 
obtain the following equation for the velocity V: 

-V+^\Vuw\^-l3,\Vuw\\ ^0 ond{uw>0}. (9) 

In one dimension, one can rigorously show that the above heuristics is valid, using the discrete- 
time (JKO) scheme introduced in ^6j (see also ^): For a given initial open interval /° and time 
step size h > 0, let us consider the sequence of open intervals /^ with i = 1,2,3, ..., iteratively 
defined by 

r+' := argmin.^op,^ interval {^rf^^^'(^^.' ^) + ^(^)} ' ^h = ^°' (10) 

where dist is a modified "distance" between open intervals (see ^ for general formulation and 
references) 

dist ((a, 6), (c, rf)) = / d{x,{a,b})dx. (11) 

J(a,c)u(b,d) 

Let Uh{-,t) and A/,.(-,t) be the associated function and lagrange multiplier to Ij^^ defined at 
discrete times t — ih. Lastly, define 

Uhi-,t),Xh{-,t) = Uh{-,ih),Xh{-,ih) ior t ^ ih < t < {i + l)h. 

Following [5|, we can then show: 

Proposition 1.1. Let u^ and I^ be as given above. 

(a) Suppose there exists a classical solution (p of 

{— A(/) < A/i (t) — 5 — (pK, cos a + XK sin a m {0 > 0} 

V <\\V(i>\'^-p-5 on^{(i>>0}r^{\D4>\^^{)} 

If h is sufficiently small depending on 5, then (p cannot cross ut from below. 

(b) Suppose there exists a classical solution (p of 

{— A0 > Xh{t) + S — (f)K cos a + xKsina in {(f> > 0} 

y > i|V0|2-/3-5 ond{(j)>0} 

If h is sufficiently small depending on 6, then (p cannot cross Uh from above. 
The proof of Proposition 1 1 . 1 1 follows that of Proposition 3.1 and Proposition 3.3 in [6 . Note 



that from ( 10 ) it follows that /^ cannot move its endpoints by more than Ch from the endpoints 
of /^~^, where C depends on the energy associated with /". It follows that along a subsequence 
Uh, \h and Ih locally uniformly converge as /i — > 0. Proposition 1 1 . 1| then ensures that the limiting 
solution satisfies (It]) and that the motion law satisfies (|9| in the viscosity sense. 

This discrete-time scheme, which relies on the gradient-flow structure of our problem, can 
thus be used to prove the existence of solution. However, it is not a very practical scheme, 



especially in one dimension. For this reason we will use a different, simpler, discrete-time scheme 
to construct classical solutions of R-Q (see Section 4.2). 



In the next section, we will state all the results proved in this paper. The rest of the paper 
(Sections Is] to [7]) will be devoted to the proof of these results. 



2 Main results and outline of the paper 

Throughout the paper, the volume of the drop Vq is fixed (we only compare drops with same vol- 
ume), and we assume that the drop is invariant in the y direction. As discuss in the introduction, 
the sliding drop problem then reduces to finding a function u(x, t) and a set W{t) — {u{t) > 0} 
such that 

u{t) = argmin{^(w) ; v G Hl{W{t)), i; > 0, v dx ^ Vo} 

r 1 1 ^ (12) 

min<^ -V + -\u^\^ - P , \u^\ } ^0 on d{u > 0}. 



2 
where V denotes the normal velocity of dW{t) and ^ is given by 



2 



/" 1 V 

■^{v) = I -\vx\ + —Kcosa — V XKsinadx. 

We note that in one dimension there is no mechanism which would allow a connected drop to 
split into several droplets (we will show this rigorously in Section 4, when we construct a solution 
using a time-discretized scheme). We thus assume that the drop is connected at time i = 0, that 
is Uiiix) satisfies {uo > 0} = (ao,6o) and 

Mo = argmin{,^(w) ; v G iIo(ao,6o), w > 0, / v dx — Vq}- (13) 



Our problem then reduces to finding W = {a{t),b{t)) satisfying (12). This leads to the following 
definition: 



Definition 2.1. Given an initial data uq satisfying ( 13 1, we say that u{x, t) is a solution of { 12) 



if there exist some Lipschitz functions a{t) and b(t) defined for t > 0, satisfying 

a(0) = ao, 5(0) = &o 

and such that 

{x;u{x,t)>0} = {a{t),b{t)) forallt>0, (14) 

Mt) 
u(-,i) =argmin{^(w); ue ffo(a(t),6(t)), u> 0, / v{x)dx^Vo}forallt>Q, (15) 

Ja{t) 

and 



min{a'(i) + i|u,(a(i))|2 - ;9(a(i)), \u^{a{t))\} = 
min{-fe'(f) + \\uMt))? - /3{b{t)), |ii.(6(t))|} = 



a.e. t > 0. 



With this definition in hand, our first result is the following: 
Theorem 2.2. Assume that x 1— > /3{x) is a Lipschitz continuous function. For any Uq satisfying 



(13), there exists a unique solution u{x,t) in the sense of Definition 2.1 



Furthermore, we have the following "comparison principle" , which plays a crucial role in our 
analysis: 

Proposition 2.3. [Comparison principle and uniqueness] Assume that x M> /3(x) is a Lipschitz 



function, and let ui , U2 be two solutions of (uM) (in the sense of Definition 2. 1 ) with support 
(ai(i),6i(t)), (a2(i),62(i)) such that 

ai(0) < 02(0) and 5i(0) < 62(0). (16) 

Then 

ai{t) < 02 (i) and bi{t) < 62 (i) for all t > 0. (17) 



Furthermore if the inequalities are strict in (I6I, then they are strict also in (171 



We stress out the fact that the solutions ui and U2 correspond to drops with the same volume 
Vo (the result would obviously not hold for two drops with different volume). The result is not 
trivial because ui and M2 do not satisfy the same equation (different Lagrange multipliers). 

Next, we investigate the existence of particular solutions characterizing the asymptotic be- 
havior of general solutions. We start with the case where (3 is constant, and we show that there 
exists a unique traveling wave solution (for a given volume and up to translation in time) , which 
is stable: 



Theorem 2.4. When /3 is constant, Problem (wM) has a traveling wave type solution of the form 
u{x,t) — v(x — ct). The profile v{x) is unique (up to translation), and the speed c is given by 

f i^o^sina ^fVo<2J^ 

c = < . (18) 



Finally, any solution of (12) converges as t —)■ 00 to a translation of this traveling wave u. 



Note that the speed c of the traveling wave is a piecewise linear function of Vb, which increases 
faster once the drop reaches a critical volume (we will see that this correspond to the volume for 
which the receding contact angle - or the gradient in the rear of the drop - becomes degenerate) . 
Note also that for small volume, the speed is independent of the relative adhesion coefficient /3. 

We emphasize the fact that when /3 is constant, one can easily show that any drop will be 
sliding down the inclined plane (no stationary solutions can exist in that case). However, the 
coefficient /? depends on the properties of the materials (solid and liquid) and it is very sensitive 
to small perturbations in the properties of the solid plane (chemical contamination or roughness) . 
In the last part of this paper, we thus consider the case of non-constant coefficient /3, and to 
simplify the analysis, we restrict ourself to periodic settings. 

Our goal is to characterize the effects of these heterogeneities on the speed of the sliding 
droplets. To this end we consider a simple case in which the constant j3 is replaced by a periodic 
function 

X H^ I3{x). 

In that case, we will see that stationary solutions do exist, at least for small volume, and 
that small enough drop may stick to the inclined plane. Large volume drops, however, will 
slide down, and the notion of traveling wave solution is replaced by that of pulsating traveling 
solutions (which are global in time solutions whose profile exhibits some periodic behavior) . More 
precisely, we will show: 



Theorem 2.5. Assume that x i~> P{x) is a periodic function. Then the following hold: 

ft) If 

max /3 — min /3 < Vqk sin a, 



then there exists a unique pulsating traveling solution u{x,t) solution of (wM) (that is a global 
in time solution satisfying u{x,t + T) — u{x + l,i) for all t). Furthermore, any solution will 
eventually start sliding with positive speed (b'{t) > for t > to) and will converge to the unique 
pulsating traveling solution as t -^ oo. 

(n) If 

max /3 — min /3 > Vqk sin a, 

then there exists S such that if the period of j3 is less than 6, then any solution will stick to the 
inclined plane. More precisely, we have 

either a{t) > a(0) - C or b{t) > &(0) - C for all t > 0. (19) 

Finally, our last result concerns the homogenization of the problem above: We now assume 
that x H-)- I3g{x) is a e-periodic function (for the sake of simplicity, we take /^^(a;) = (3{x/e)). 
First, we state the following lemma which will be proved in the appendix: 

Lemma 2.6. Assume that x i— >■ P{x) is periodic. Then, for all q CzR, the equation 

x'{t) = q-p{x{t)) 
has a unique global solution (up to translation in time), which is periodic in time. We define by 

r{q) := lim ^ ' ^ ' (20) 

t— >-oo t 

its effective speed. The function q^-^ r{q) is monotone increasing and continuous, and satisfies 

r{q) = 0, for q G [min /3, max /3]. 

Note that the fact that r(g) can be zero in a non trivial interval is a classical and important 
aspect of the homogenization of contact angle dynamics (see Figure [I]) . It typically implies 
pinning of the contact line: for related results in higher dimensions, we refer to [5j and [TU] . 

Wc can now state our homogenization result: 



Theorem 2.7. Let u'^ be a solution of (12) with initial value uq and where /3 — /3{x/e) with (3 a 



periodic function. When e goes to zero, u'^ converges uniformly to a function u{x,t) solution of 

a.e. t> 0. (21) 



(14|-(15l satisfying the following velocity law: 

r min{a'(t)+r(i|u,(a(i))|2),|u,(a(t))|} = 
1 min{-&'(t) + r{}\uMtW), l«x(&W)l} = 



Furthermore, there is at most one solution u of the homogenized problem ( 14)-( 15 )-(21 1 with 
given initial data uq. 



To conclude this analysis, we note that the homogenized problem described in Theorem |2.7| 
has traveling wave solutions for all volume (for small volume, the traveling wave has speed zero, 
so it is really a stationary solution). The speed of that traveling wave can be determined as 
follows: 




Figure 1: Homogenized velocity function r(q) as a function of the slope q when (3 is given by 
/3(x) = 1 + 0.3 sin a;. The straight line is the curve y = q — {/3) 

• If r(VoKsina) + r(0) < 0, then there exists qo such that 

''(90 + VqK sin a) + r{qn) = 
and the speed of the traveling wave solution is given by 

c = r{qn + VoKsina) == -r{qn) 
(note that go may not be unique, but c is unique). 

• If r{VoKsma) + r(0) > 0, then the speed of the traveling wave solution is given by 

c = r(VbK sin a) 

(when this happens, the gradient in the rear of the drop of the traveling wave solution is 
zero) 



Of course, we recover the speed of the traveling wave for constant coefficient /? (formula (18l) 
when we take r{q) — q — /3. Figure ^ compares the speed of the traveling wave as a function of 
VoKsina for the homogenized problem ( |21| when /3g(a;) = 1 + 0.3sin(a;/£) (i.e. with the function 
r shown in Figure ll} and for the constant case (3 — {j3^) — 1. Note the Lipschitz corner near 
VoKsina = 2 which corresponds to the critical volume for which the gradient in the rear of the 
drop becomes zero. 

The rest of the paper is devoted to the proof of these results. 



3 Preliminaries 



3.1 Another formulation for (12) 



Definition 2.1 gives the most natural notion of solution for (12). However, it is not the most 



practical one. In this section, we derive an equivalent formulation which will be more convenient. 




0.5 1 1.5 2 2.5 3 3.5 4 



Figure 2: Speed of the traveling wave solution as a function of VoKsina for the homogenized 



problem (21) when I3^{x) = 1 + 0.3sin(a;/e) (solid line), and for the constant case /3 = {Pe) = 1- 



First, we note that if u is a solution in the sense of Definition 2.1 then u{-,t) solves: 

—Uxx — M^) ~ UK cos a + {x — b{t))K sin a in (a(i), b{t)) 

u{a{t),t) ^u{b{t),t) =0, 

udx = Vq. 



(22) 



It is tempting to replace (14)-p5| by (22); However, solutions of (22) might take negative 



values (for large intervals (a(i), 6(t))). We thus write: 



Definition 3.1. Given an initial data Uq satisfying ( 13 ), we say that u{x, t) is a solution of { 12) 
if there exist some Lipschitz functions a{t) and b{t) satisfying 

a(0) = ao, 5(0) = bo 

such that u{-,t) solves (pM) for all t > and a{t), b{t) satisfy 

r mm{a'{t) + i|«,(a(i))|2 - /3(a(i)), u,(a(t))} = 

1 b'it)='^\uxibm'-Pibit)). 



(23) 



In this second formulation, the positivity of u in (a(t),6(f)) is enforced by the condition 
Uxici) ^ 0- Indeed, we can show: 



Lemma 3.2. Let u{-,t) be the solution of [22), then —Ux{b{t),t) > 0. Furthermore, u{x,t) > 



in {a(t),b{t)) if and only Ux{a{t),t) > 0. In particular, Definitions 2.1 and 3.1 are equivalent. 



Note that it might have seemed more natural to write the second equation in ( 23 ) as 



min{-6'(i) + \\ux{b{t))\^ - P{b{t)),-Ux{b(t))}- 



However, lemma 3.2 shows that the condition —Ux{b{t)) > always holds. 
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Proof of Lemma \3.S\ We drop the t dependence in what foUows. Smce Vb > 0, there exists 
Xq e {a,b) such that u{xq) = max-x(z(a,b) u > 0. Since —u"(xq) > 0, (22) imphes 



A + {xq — 6)k sin a > u(xq)kcos a > 0. (24) 



Assume now that u has a local minimum at a point xi G (a, b). Then —u"{xi) < 0, and so (22) 
implies 

A + {xi — 6)Ksina < u(a:i)KCOSQ!. 

Since u(xi) < v{xq), we deduce 

Xl < Xq, 

and SO we always have u{x) > in {xQ,b). In particular, we have —u'{b) > 0. In order to show 



that —u'{b) > 0, we note that (24 1 also implies A > 0, and so the function h{x) = —fJ.{x — b) 
satisfies 

^h" < X — hn cos a + {x — 6)k sin a = A — [fin cos a + k sin a] (b ~ x) 

in a neighborhood of b. A Hopf's Lemma - type argument yields —u'{b) > 0. 

Assume now that u'{a) > 0, and that there exists Xq € (a, 6) such that u{xo) < 0. Using the 
boundary conditions, we deduce that there exist at least three distinct points in [a, b] where u' 
vanishes, and so there exists at least two distinct points j/i, y2 in (a, 6) where u" vanishes. 

Therefore, the function w = uk cos a — X — {x ~ b)Ksma, which solves 

w" — Kcosaw in (a, 6) 

is such that w"{yi) — w{yi) = 0, i — 1, 2 and the maximum principle implies w = in (2/1,2/2) 
and thus w = in (a, &), which is impossible (look at the values of w at a and b). 

We have thus shown that u'{a) > implies that u > in (a, b). Conversely, if u > in (a, b), 
it is readily seen that u'{a) > 0. 

□ 

Finally, we give a third definition which will be very useful later on. We will show in Propo- 
sition [3^ (iw) that there exists a critical length £c > (depending on the volume Vq and the 



inclination a) such that the solution of ( 22 ) is positive in (a(t), b(t)) if and only if b{t) — a{t) < ic- 
This length f^ is the longest possible length of the support of the drop. This leads to the following 
third formulation: 



Definition 3.3. Given an initial data uq satisfying ( 13 1, we say that u{x, t) is a solution of { 12) 



if there exist some Lipschitz functions a{t) and b(t) such that u{-,t) solves {22) and 

a(0) = ao, KO) = &o 

with 

r min{a'(t) + i|u,(a(i))|2 - /3(a(i)), a{t) - b{t) + 4} = 

\ b'{t)^\\uMt))?-p{m), 

a.e. t > 0. 

Equivalently, denoting £{t) = bit) — a(t), we can write: 

r max{f (t) - [i|u,(6(t))p + i|u.(a(i))|2 - /3(a(i)) - mM ,m - 4} = 

{ b'it)^l\uMtW~P{b{t)), 



(25) 



(26) 



where Ux{b{t)) and Uj.{a{t)) only depends on l(t) (see Proposition 3.5 {Hi)). 
We then have: 
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Proposition 3.4. Definition \3.1\ and \3.3\ are equivalent. 

In order to prove this proposition, we will need to derive several important properties of the 



solution of the obstacle problem (15), which is done in the next section. Proposition 3.4 will be 



an immediate consequence of Proposition |3.8| below 



3.2 Analysis of the obstacle problem (15) 



Throughout this section, we fix a volume Vq, and all constants will depends on Vq, a and k. 
Given a < b, we consider u solution of 

u = argmin{^(u) ; v E H^ia, b), v>0, v{x) dx = Vq}. (27) 

J a 

Then classical arguments for obstacle type problems yield that u E C^'^{a, b) solves 

min{— m" — A — (x — 6)Ksina + MKC0sa, u} = in (a, 6), u(a) = u(fe) = 0, udx = Vo (28) 

(where A G M is determined by the volume constraint). 

The main properties of u are gathered in the following proposition: 



Proposition 3.5. Let u be a C^'^ solution of (28), th 



(i) {u > 0} is an interval of the form (a', b) and —u'{b) > 0. 



(ii) u is the unique solution of \2^ , and is also solution of (21). 

(Hi) X, — |{m > 0}|, ~u' [b) and u'{a') only depends on a and b through (b~ a) and are decreasing 
functions of b — a. 

(iv) There exists ^^ > such that {w > 0} = (a, h) if and only ifb—a < £c. Furthermore, u'{a) > 
if and only if b ~ a < £c- 

The critical length (^ is the longest possible length of the support of the drop. We recall that 
it depends on Vq, k and a. It will play a very important role in the analysis of our problem. 



In view of (Hi), we can define 



I '/ m2 rj/i, \ 1| I n,\\2 



G{b-a):^-\u'{a)\\ H {b - a) := -\u' {b^ (29) 



where u is the solution of (28). The functions £ i— > G{£) and £ i~^ H{£) are then both monotone 



decreasing, and multiplying the equation (28) by u' and integrating over (a', 6) yields 



^u'{b)^ - ^u'ia'f ^ KsinaVo 

or 

H{£)-G{£)^KsmaVo. (30) 

We also define 

F:£^G{£) + H{£), (31) 

which is a decreasing function as well. 

Finally, we have the following lemmas, which will be useful later on: 
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Lemma 3.6 (Lipschitz regularity). For all Sq > 0, there exists C > such that 

< H{eo) - H{io + Tl)<Crj 

and 

0<G{eo)-Gi£o + v)<Cv 

for £q > 6q and rj > 0. In particular, it follows that H and G (and therefore F) are locally 
Lipschitz on (0, oo). 

Lemma 3.7 (Non-degeneracy). For all (5o > there exists c > such that 

H{eo) - H{io + v)>crj (32) 

for all £o < ic ~ ^a ind rj such that io + rj < ic 

Next, we consider the following boundary value problem with volume constraint: 

{—v" = X — VK cos a + (x — b)K sin a in (a, b) 
(33) 
v{a) = v{b) = and J v dx = Vq. 

where A is a Lagrange multiplier for the volume constraint. Then, we have: 



Proposition 3.8. Let v be the solution of (33). Then there exists (.c depending on Vq, a and k 
such that the fallowings are equivalent 

(i) v'{a) > 



(ii) V > is also the solution of (27) 
(a) b — a < £c 
Furthermore, v'{a) > if and only if b — a < £c. 

This proposition implies Proposition |3.4[ 

The remainder of this section is devoted to the proof of these results. 



Proof of Proposition \3^ 
(i) Let (c, d) be a connected component of {u > 0}, then 



— u" = X + [x — b)K sin a — uk cos a (34) 



in (c, d). Multiplying (34 1 by u' and integrating over (c, d) yields: 



u'{d)=u'{c)+ 2k sin a udx>0. 

li d y^ b, then the regularity of u implies u'idY = 0, a contradiction. So d = 6, |m'(&)| > and 
(c, 6) is the only connected component of {u > 0}. We also have the following formula: 

u'{bf -u'{a'f = 2KS\naVQ (35) 
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(ii) To prove the uniqueness of the solution of (|28[), we assume that ui, U2 are two solutions of 



(281 (with Ai and A2 respectively). We know that {u^ > 0} = (a^, b) and we can assume (without 
loss of generality) that a < a2 < a'^ . 
li V — ui — U2, then 

~v" — Xi — \2 — VK cos a in (a'^ , b) . 

Furthermore, v{a'i) < 0, v{b) — and J , v dx > 0. So u reaches its positive maximum value 
somewhere in {a'i,b). We deduce 

Ai > A2. 

Next, we note that since 

—Ui > Ai + (x — 6)Ksina — uiKCOsa in {aj, 6}, 

we also have 

—v" > Ai — A2 — VKCOsa > —vKCOsa in (a'j, 6), 

with ^(02) = v{b) = and J , v dx = 0. If w ^ 0, then v has a (strictly) negative minimum value 
in (02,6), a contradiction. We deduce v = and so ui — U2- 



(iii) Equation ( 28 ) is invariant by translation, and thus the uniqueness obtained in (ii) implies 
that u{- -\-b) and A only depends on £ = ib — a) (in particular, u' [a) and —u'{b) only depends on 

(■)■ 

The proof of (ii) above also implies that A is a decreasing function of [b — a). We now take 



a2 < ai < b and let Ui (resp. U2) be the solution of (28) on the interval (ai,&) (resp. (02,6)). 
The function v = ui — U2 satisfies 

v{ai) < 0, v{b) = 0, and v>0 

SO there exists a point xq G [ai, b) such that w(a;o) — 0. Since 

— w" + WKCOsa = Ai — A2 > in(a;o,6), 
the maximum principle gives w > in [xq, b). In particular, —v'{b) > and so 

-U[{b)>-U'2{b). 



Finally, if ^2(^2) = 0, then W2('^2) < '"U'^i)' while if ^2(^2) > 0, then (351 implies 

^2(02)^ = ""'2(6)^ — 2KsinaVb 

and 

u[{a[f = u[{b)^ - 2Ksinayo > u'2ibf - 2Ksinayo = '«2(a2)^ 

hence a[ = ai and ui(ai) > M2(a2)- 



(iv) Note that if {u > 0} 7^ (a, 6), then u is the unique solution of (28) in any interval larger than 
(a, b). So we can define 



= min{£ ; the solution of (28) in (0,^) satisfies {u > 0} ^ (OJ)} 
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We only need to show that this set is non empty (that is ic < c»), and that this is equivalent to 
u'{a) 7^ 0. 

Assume that {w > 0} 7^ (a, h). Then 

— u" < A + (x — 6)k sin a in (a, b) 

and so 

—u'{h) + M'(a) < A(6 — a) (6 — a)^KsinQ; 

Since u'{a) > and u'(6) < 0, we deduce 

{b — a)Ksina < A 

Since A is a decreasing function of (6 — a), this implies that [b — a) < C for some constant C . 

Finally, we note that if fe — a < ^c, then u' {a) > 0, while the regularity of u implies u' {a) — 
when b — a > Ic- 

□ 

We now turn to the proof of the Lipschitz regularity of the functions G and H: 
Proof of LeTnTna \37S[ Note that H and G are constant for £ > £c, so we can always assume that 

So<eo<£o + ri< 4. (36) 



We now consider the solution u of (28) with a = and b = £q and we denote by A(^o) the 
Lagrange multiplier appearing in (28). We recall that A is a decreasing function of ig and we 
claim that 

Xiio)-Xi£o + v)<Cv (37) 

for some constant G depending on Sq. To prove this, we consider v solution of the same equation 
as u (with the same A) but on the interval (0,^o + v) instead of (0,£o): 

—v" = X{£o) + {x — £o)Ksina — vKCOsa in(0,£o + '7), v{0) ~ v{£o + -q) — 0. 

We note that v > u and that the function h = v ~ u solves 

hn cos a — /i" = on (0, ^o)- 

In particular, on the interval (0,-^o)j h must take its positive maximum at a; = £q. Now on 
{£o,£o + rj) we have h ^ v < Giij for some Ci > with Ci depending on Sq- We conclude that 
h < CiV in (0, £0 + v) ^^d thus 

"" V dx < Vq + C[r] (38) 



for some constant C( depending on Sq and £c (recall that £0 + 1] < ic)- 



Next, we consider w the solution of (28) with a = and b = io + rj (and so A = X{£o + rj))- 
The function g = v — w solves 

gK cos a — g —a 

in (0,^0 + v)7 "^ith 5(0) ~ g{£o + v)=0, where (since A is decreasing) 

a = X{£o) — X{£q + rj) + rjKsina > 0. 
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Finally, we note that the solution of — (/>" = /ia in (0,^o + v) with zero boundary condition 
satisfies <^KC0sa — 0" < cr, if /i is small enough. The comparison principle thus implies that 
g > (j) and so 

gdx — I V dx — Vo > I (pdx > Cfia 



for some small /i > 0. Equation (38) implies 

C[r]> Cfia 

and so 

A(£o) - A(4 + v)<Cr, 



which proves (371 



We can now prove the Lipschitz bounds for the functions H and G by comparing u (solution 
of (28) on (0,£o)) and w (solution of (p8| on {0,£Q + r])): The function (/(x) ^ u{x) — 'w{x) solves 



qK cos a — q" = A(£o) — X{£q + rj) + r/K sin a, in (0, ig) 
with q{0) = and g(4) = -w{£q) < 0. So (l37| implies 



qncosa-q" <Ct] in (0,4), with g(0) = and g(£o) < 0. 

We deduce q'{0) < Crj and so 

u'(0) < w'{Q) + Crj 

which proves our statement for G. The proof for H is similar (fixing 6 = and choosing a — —£o 
or a = —£q — rj). D 



The proof of Lemma 3.7 follows from similar argument: 



Proof of Lemma 3. ? We recall that £q < £c — 5q and we fix fe = 0, 02 = — (^0 + v) and ai ~ —£q. 
We note that it is enough to prove (32 1 for 77 < {£c — 4)/2, so that we can assume that to+rj < 
4 - <5o/2. 



We now consider ui (resp. U2) the solution of (28) on the interval (ai,0) (resp. (02,0)). The 
function h = ui — U2 satisfies in particular 

h{ai) < 0, h{0) = 

and we are going to show that max(aj 0) h> cq for some c > 0. 
First of all, the volume constraint implies that 



ph pa 

/ hdx = 



U2dx > 



(39) 



Furthermore, since Iq + rj < £c — 6q/2 the monotonicity of G implies 

u'^ia2)>^2G{i,-do/2)>0 
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and so ^2(01) > crj for some c > 0. Also, Lemma 3.6 gives 



/i'(ai) — u'i{ai) — u'2{ai) 
= G{eo) ~ u'^iai) 

< G{£o + r^) + Cv~u'2ia,) 

< C?7 + ^2(02) - U2("i) 

< Ct] 

We deduce that there exists c > such that 

h{ai) — ~U2{ai) < —crj and /i'(ai) < crj. 



Finally, we have —h" + hK cos a = Ai — A2 in (ai, 0), and so (37) implies 

< -h" + hncosa < C-q. 



(40) 



(41) 



Consider now g solution of —g" + gncosa = C on (ai,0) with g{ai) = — c and g'{ai) = c. 



Then (40) and (41 1 imply that h < ijg in (ai,0). 

Let now yg € (ai,0) be such that 5 < in (ai,2/o) and g{yo) = (yo only depends on the 
constant c and C appearing in (40) and (41|), then 



yo ryo 

hdx < 7] gdx < —cq 



(42) 



with c > 0. 



Equations ( 39 1 and ( 42 ) imply that 



h{x)dx > cq 



vo 



and so there exists a point Zq G {yo, 0) where h{zQ) > cq. A Hopf's lemma type argument (using 
(41)) now yields |/i'(0)| > cq and so 



l«'i(0)|> 14(0)1 +C77. 



The result follows. 



□ 



Proof of Proposition 3. 



(i)=4>(ii) If v'{a) > 0, then Lemma 3.2 implies that w > in (a, b), and so v is also the solution of 
p8|, and thus of &7\. 



(ii)=^(iii) If V is solution of (33 1, then w > in (a, b), and since v cannot be zero a an interval of 
(a, b), we must have {v > 0} = (a, b). Proposition 3.5 (iv) gives b — a < £c- 

(iii)=>(i) Finally, ii b — a < £c, then the solution u of (p7| satisfies {u > 0} = (a, b), and is thus 



also a solution of (33). If follows that v'{a) — u'{a) > 0. 



□ 
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4 Proof of Theorem 12.2 



4.1 Uniqueness and Comparison principle 

In this section, we prove the comparison principle, Proposition [231 which also implies the unique- 
ness part of Theorem |2.2[ 



Proof of Propositi on \2. 3| Thanks to Propositio n |3.4[ we can assume that u is a solution in the 
sense of Definition 3.3 Using the notations of (29), we can rewrite (25) as follows 



min{a'(i) + G{(.{t)) - l3{a{t)),£c - iit)} 
h'{t)=H{t{t))~P{b{t)), 







(43) 



with (. I-)- -ff (£) and £ i->- G{t) monotone decreasing. We now define 

fit) = max{ai(t) - a2{t)M{t) " &2(i)}- 

Since the functions ai{t) and bi{t) are in PF;q'^(0,oo), f{t) is in W^;J^(0,oo) and we have 

r a'l {t) - 4 (i) if ai{t)- 02 (t) > bi {t) - ^2 (t) 

f'it)=l b[{t)-b'2{t) iia,{t)~a2{t)<biit)-b2{t) 

[ b[ (t) - b'2 (t) = a[ (t) - a'2 (t) if ai (t) - 02 (t) = &i {t) - 62 (t) 

and so 

r a[it) - a'2{t) a i2it) > hit) 

1 b[{t)-b',(t) if4(t)<£i(i) 



a.e. t > 



When £2{t) > ii{t), we have in particular £i{t) < £c, and so (43) gives 



We deduce 



fit) < 



a[(t) = -G{£i{t)) + /3(ai(t)), 4(t) > -G{£2{t)) + I3{a2{t)). 
" G(4 W) - G(Mt)) + /3(ai(i)) - I3{a2{t)) if £2 W > hit) 



Hi£,{t)) - H{£2{t)) - P(b^{t)) + P{b2{t)) if £2(0 < h{t) 

and the monotonicity of G and _ff implies 

r if |ai(i) - a2(i)| if ai(i) - aalO > 61 W - b2{t) 
^^' - \ K\bi{t) ~ b2{t)\ if ai(t) - a2{t) < biit) - b2{t) 

with K = ||/3'||loo. We thus have 

f'{t)<K\f{t)\ a.e. t>0. 

In particular, /^(t) < X/+(t) and so f+{t) < /+(0)e^*, so if /(O) < 0, then f{t) < a.e. t > 0. 

We also have /i(i) > -Kf'_{t) and so /_(i) > /_(0)e--^*, so if /(O) < 0, then f{t) < a.e. 
t > 0. The resuh follows. D 



18 



4.2 Existence of a solution 



Observe that, when /3 does not depend on x, formulation (26) imphes that the problem can be 
reduced to solving first an equation for £{t), of the form 



max{/(t) - [F{e{t)) - 2/3],£{t) - 4} - 0} 



with F defined by (31 1, and then an equation for b{t) of the form 

b'{t) = H{e{t)) - p. 



(44) 



(45) 



When /3 depends on x however, one cannot decouple the equations for a and h. We will thus use 
a discrete scheme to prove the existence of solutions in this general framework. 



4.2.1 Discrete-time scheme 

To prove the existence of a solution for general coefficient /3, we introduce a simple discrete time 
scheme (simpler than the gradient flow scheme described in the introduction): For h > {) small, 
we describe the evolution of the drop as follows: 

Assume that the support of the drop at time t„ = nh is the interval (a„, 6„) and denote by 
Un{x) the corresponding solution of (27 1. We define 



"n+l/2 — "r: 
^n+1/2 = bn 



I3{ar. 



K(an)P 



i|<(6„)p-/3(6„) 



h 



Next, we define u^+i as the solution of the obstacle problem (27 1 with a = an+i/2 and b = 6,1+1/2- 
Proposition 3.5 (i) implies that there exists a„+i and fe„+i = 6n+i/2 such that 



{un+i > 0} = (a„+i,6„+i). 

This scheme defines a sequence of functions {u„(a::)}„gN and a sequence of intervals {(a„, 6„)}nGN- 
We then define u'^{x,t), a^{t) and b^{t) continuous piecewise linear functions such that 

u '(x, nh) := Un{x) 

a^{nh) :— a„, b'^[nh) :— bn- 
Our goal is now to pass to the limit /i — >■ 0. 

First, we observe that since fe„+i = 6„+i/2, we have 

bn+l -bn I. 



K{bnr-P{bn). 



(46) 



On the other hand, we only have a„+i > a„+i/2 and if a„_|_i > a„+i/2, then u'„^]^(a„-|_i) = 0. So 
we can write 



, > /3(a„) - -|uUan)P, with equality if u'„+i(a„+i) > 0, 



which implies 



mm < "'"^\ + 2 K(a«)P - /3(a„),<+i(a„+i) 



h 



Finally, using Proposition 3.8 we can also write: 



O-n+l 0-n t I , N|2 of \ 

mm< ^ ^\Un[an)\ -p{an),a„ 



(47) 



(48) 
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4.2.2 Limit h^O 

We now study the limit /i — > of the discrete model introduced above. Note that u„ satisfies, 
for some A„, 

—u'^ — Xn ~ UnKCOsa + {x — &„)Ksina in (a„, 6„), 

Unian) = Un{bn) = 0, (49) 

, JUndx = Vq. 

We then have the following proposition: 

Proposition 4.1. For all n, the following holds: 

(i) |u^(a„)| and |u^(&„)| o.^^ uniformly hounded. More precisely, there exists a constant M 
depending only on k, a, Vq and the initial data such that 

|<(a„)|<|<(6„)|<Af. 

(ii) There exist constants £ > and A depending only on k, a, Vq such that bn — a„ > | and 

|A„|<A. 

(Hi) a^{t) and b'^{t) are uniformly Lipschitz continuous with respect to h and t. 

Before proving this proposition, we show that it implies the existence of a solution (Theo- 



2.2) 



Proof of Theorem\2^ Proposition 4.1 implies that a'^{t), b'^{t) and £''(i) = b'^{t) - a'^{t) are 
Lipschitz continuous uniformly with respect to h and that there exists | independent of h such 
that 

e<i''it)<ic Vi>0. (50) 

In particular there exists a subsequence /i — >■ and some Lipschitz continuous functions a(t) and 
b{t) such that 

a'^{t) -^ a{t), b'\t) -^ b{t) uniformly with respect to i G [a, b] 

for all a < 6 e M. 

Next, we note that for any i > 0, u^{x, t) solves 

~Uxx{-,t) = A(^''(i)) - KCOsau+ {x - 6''(t))Ksina, in (&'',&'*) 

and so the function v'^{x,t) — u^{a^{t) + x£'^{t),t) solves 

-Vxx{-,t) ^t{tf [\{t{t)) - Kcosav + t{t){x ~ l)KSina\ , in (0,1) 

with 

i;''(0,i)=w''(l,i)=0. 

In particular, it is readily seen that x i~^ v^{x^ t) is bounded in C^(0, 1) uniformly with respect 
to t, and that v^ and w^ are Lipschitz continuous with respect to t, uniformly with respect to h. 
Finally, v^ converges locally uniformly (with respect to x and t) to a function v solution of 

-v.^^{-^t) ^ i{tf [\{£{t)) - KCOsav + e{t){x- l)Ksina] , in (0,1) with i;(0,i) = v{l,t) =0. 
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Writing 



Am)^.'(^,. 



we deduce that u^ converges locally uniformly (with respect to x and t) to a function u solution 
of 

—Uxx{', t) — X{i{t)) — K COS au + {x — b{t))K sin a, in (a(t), b(t)) 

satisfying 

u{a{t),t) ^u{b{t),t) = and u{-,t)dx = Vo. 



Furthermore, u^{a {t),t) and u^ib {t),t) are Lipschitz continuous function in t, uniformly 



with respect to h. We can thus pass to the limit in (46 1 and (48). For instance, we note that 



(46) implies 



for t e (nh, {n + l)h), and so the Lipschitz continuity in time implies 

b''\t)^^\u'i{b\t),t)\^^(3{bHt)) + O{h)a.e.t>0 



Passing to the limit < — > and repeating this argument with MSl yields ( 25 ) 



□ 



Proof of Proposition 4-1 For the sake of clarity, we drop the index n when no ambiguity is 
possible. 



(i) In view of Poposition 3.5 we only need to show that the length of the support b ~ a cannot 



be too small: First, if u'(o) = 0, then Proposition 3.5 (iv) implies that b — a = ic, and (30) gives 

-|u'(6)| =VbKsina. 



Next, if u'(a) ^ 0, then (46) and (47) implies that the length of the interval 6— a at the next 
time step will be given by 



b- a + 



u'ibf u'{a) 



'^„^2 



/3(a„) - fi{bn) 



> b-a 



i(U\'i 



u'{b) 



P{an) - /3(6„) 



In particular, if |w'(&)| > 2sup-\//3, then the length of the interval increases and therefore |w'(6)| 
decreases. We deduce that if u'{a) ^ 0, then 



and the result follows. 



W{a)\ < |«'(6)| <max{2supV/9,|u^(6)|} 



(ii) The proof of (i) above clearly implies that 6 — a is bounded below. Furthermore, integrating 
the equation satisfied by u over (a, 6), we get: 



X{b-a) 



— ?i'(a) — u'(6) + Vqk sin a {b — a)Ksina 

< 2M + VbKsina. 



The result follows. 
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(iii) The discrete motion law ( 46 ) and (i) imply that there exists a constant C > such that 

^"+^ ^ ^" <C foralln>0. 

h 



Similarly, if |M^_|_;^(a„+i)| ^ then (47) and (i) imply 

< C for all n > 0. 



O-n+l ~ 0,n 



h 

Finally, if u^_,_]^(a„+i) — 0, then bn+i — fln+i = ^c and since bn — an < ic we have 

fln+l ^ In bn+1 — b„ 



< 



<c. 



and ( 47 1 implies 



h 

a-n+l ~ a„ 1 2 ^ 1 ,^2 



>/3(a„)~-M^>--M^ 



The result follows. 



□ 



5 Asymptotic behavior when (5 is constant: Existence of 
traveling wave solutions 

In this section, we study the long time behavior of the drops when f3 is constant (homogeneous 



media), and prove Theorem 2.4 



As mentioned earlier. Problem ( 12 ) is much simpler in this case since the length of the support 
e(t) = b{t) - a{t) solves 



max{£'(t) - [F{(.{t)) - 2(3], £{t) - 4} = 



(51) 



with F defined by (31 ). Once the solution of (51 1 is found, the solution of ( 12 ) is fully determined 
by solving 

b'it) = Him)-/3. 



A traveling waves type solution of (p2|-(25) is a solution of the form 

u{x, t) — v{x — ct). 
For such a solution, the endpoints of the support (a(t), b{t)) of u satisfy 

a' {t) ~ c and b' {t) — c for all t. 



(52) 



and the length b{t) — a{t) is constant equal to some Iq. Equation (51 ) thus reduces to finding £q 
such that 

max {- [F{£o) ~ 2/3] , £„ - 4} - 0. (53) 

To solve this equation, we recall that F is continuous and monotone decreasing on the interval 
(0,4]- So we have the following: 

1. Either F{£c) < 2/3, in which case, there exists a unique £q < 4 such that F{£q) = 2/3. 

2. Or F(£c) > 2/3, in which case we have £q — £c- 
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This proves the existence of a unique travehng wave, and we now derive the formula ( |18[ ) for 
the speed c: Let a and h be such that b — a = io and let u be the corresponding solution of ([T]) 



with support (a, &). li Ux{a,) > , then (25) and (52) imply 



and so ( 30 ) yields 



This is possible only if /? - 



^\ux{a)\'^ > 0, that is if iVoKsina < /3. 



liu^ia) = 0, then (52), (25) and (30) imply 



c ~ H{£q) — (i ~ VoKsina — /3. 



This implies formula (18). 



Finally, the stability of £o is an immediate consequence of the monotonicity of F since it 



implies that any solution of (51) satisfies 



f(t) <Oif £<4 
i'{t) >Oif £>4- 

and so £(t) — >■ £o as i — J- oo. Since the profile u{x,t) is completely determined by the length of 
the support, the convergence of £{t) implies the uniform convergence of u(x, t) to the profile of 
the unique traveling wave. 



6 Asymptotic behavior when P is periodic: Existence of 
pulsating travehng solutions 

In this section, we investigate the long time behavior of the solution when the function /? is 



periodic, and prove Theorem |2.5 



This is of course more delicate than the case where B is 



constant, since we cannot reduce (12) to a single equation for the length £{t) of the support. 



First, we observe that for any solution we have (using (30)): 



a'(i) + 6'(i) 



> 



> 



iKibm'-lwAam' 

Vqk sma + /3(a(t)) - /3(6(t)). 



/3{a{t)) - mt)) 



In particular the condition 



max /3 — min /3 < Vqk sin a + S for some 6 > 



(54) 



(55) 



guarantees that a'{t) + b'{t) > 6 for all time. This implies that no stationary solution can exist 
under condition (55), and that any drops must slide down the inclined plane. 



In order to prove the first part of Theorem 2.5 (existence and stability of Pulsating traveling 



solutions), we will prove that (55) actually implies that b'(t) > for all (large enough) time, 
and that we can thus rewrite the equations for a{t) and b{t) using x = b{t) instead of t as a 
parameter. We thus start with the following lemma: 
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Lemma 6.1. The function 1 1— > b{t) is in C^'^ and if (55) holds then there exists r] > such that 
b'{t) > T] at least for t > Tq ^ 2^. 



Proof. Let 77 G (0,(5/8) be a small number, to be chosen later. Equation (54) and condition (55) 
implies 

a'{t) + b'{t) > 6, for all t > 0. (56) 

So as long as b'(t) < rj, we have a'{t) > 5 — rj and thus 

£'(t) = b'{t) - a'{t) <-5 + 2ri< -(5/2. 

In particular, we cannot have b' [t) < 77 in (0,ro) (with Tq — 2jYm), or the length of the droplet 
would shrink to a point. 

Therefore if the lemma is false, then there exists tg > ti such that &'(to) = "H ^-nd b'{t) > rj 
for t G {ti,to). Next, we note that t 1— ?> £{t) and t H> b{t) are Lipschitz continuous and so (43) 
implies that t 1— >■ b'(t) is also Lipschitz continuous. Therefore, there exists h < ^ (where K is 
the Lipschitz constant of b') such that 

b'{t)<2ri forte (to- /i, to). (57) 

Proceeding as before (using the fact that 77 < (5/8), we deduce that 

e'{t) < -5/2 forte (to -ft, to) 
and so 

e{to) < £{tQ -h)--h 

(in particular, t{t) < £c and so Ux{a{t),t) > for t e (to — h,to)). Lemma 3.7 implies 

H{e{ta)) > H(i{to - h)) + coSh (58) 



for some cq > 0. Finally, (57) implies 

\b{to) - b{to - h)\ < 2hr] 
and so ([43| gives 



77 = 6'(to) = i/(4) - /3(KM) > Hie{ta-h))+co6h~l3{b{to~h)-\\l3'\\oo2hT] 

> b'{ta-h)+ca6h-\\p'\\^2hr, 

Since we can always choose h such that 6' (to — h) > rj (by taking h < to — ti), we obtain a 
contradiction if we choose 77 such that 



V< 



Co (5 



□ 



Proof of Theorem 2.5 (i). Lemma [6.1| iinplies that 1 1— > 6(t) is bijective from (To, 00) onto (6(To), cxd). 
We can thus re-parametrize the motion of the drop using x = b{t) as our new parameter. We 
now write the equation satisfied by y(x) = £{b~^{x): 

max{y\x)b\b-\x)) - [F{y{x)) - Pia{b-\x))) - f3{x)],y{x) - 4} = 0. 
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Using the fact that b'{b ^{x)) > and a{t) = b{t) — £{t), the equation can be written as 

niax{y'(x) - ^7(^^r^ [F{y{x)) - I3{x - y{x)) - P{x)U{x) - 4} = 0. 

Finally, using the equation for b'{t), we deduce that y solves 

max{y'(x) - F{y{x),x), y{x) - 4} = 0. (59) 

with 

X-. F{y{x))-P{x-y{x))^p{x) 

^'^' H{y{x))-P{x) 

In particular, the function (a;, y) i— )■ J-{x, y) is periodic with respect to x. 

We now need the following lemma which will be proved later on: 
Lemma 6.2. There exists a 1 -periodic function z{x) such that 

lim y{n + x) — z{x). 

In order to conclude, we rewrite the equation for b(t) as 

b'it) = Hiyibit)))-/3{b{t)) 
= H{b{t)) + ip{b{t)) 

with T-L{b) := H{z{b)) — (3{b) periodic function, and Lp{b) := H{y{b)) — H{z{b)) satisfying 

199(6)1 — > as 6 -^- 00. 

Since b'{t) > rj, there exists a sequence t„ — >■ 00 such that 6(i„) = n. Let us consider 
bn{t) '■= b{tn + t) — n, which solves the following equation: 

b'^it) = n{bn{t)) + ifibnit) + n), bniO) = 0. 

We deduce 

\bn+k{t) - bnit)\ < e^*(|v?(n)| + |^(n + fc)|), 

and thus {&«} is a Cauchy sequence: it converges to a monotone increasing function bao{t). Since 

bn{tn+l - tn) = ^n+l(0) + 1, 

we deduce that i„+i — i„ converges to the unique to satisfying 600 (io) = 1 and that feoo is 
io-periodic. 

Hence, we can show that b{t) — b^{t) converges to as t goes to 00, and defining foo(i) = 
z{bao(t)), and aoo(t) = booit) — £00 (t), we obtain a pulsating solution {aooit),boo{t)). 

Lastly, the uniqueness of the pulsating solution with positive speed (which is guaranteed by 



Lemma 6.1) follows from the comparison principle (Proposition 2.3) 



□ 



25 



Proof of Lemm,a \6.^ We recall that y{x) is a bounded solution of (59). For ?i e N, we denote 



yn{x) = y{x-\-n). The periodicity of J^ with respect to x implies that it is also a bounded solution 



of (59) with initial condition y„ := y{n). 

If yi = yoi then the uniqueness implies that y(a:; + 1) = y{x) for all x and we can take 
z{x) = y{x). Otherwise, if, for instance yi > yo, then the uniqueness implies that yi{x) > yo(x) 
for all a; > and so y2 = yi(l) > J/o(l) ~ yi- We deduce y2{x) > yi{x), and iterating this 
argument, we get 

yn+i{x) > yn{x) for aU n g N. 

So {yn{x)} is monotone and bounded, and thus there exists z{x) such that lim„_j.oo yn{x) — z{x). 
Since j/„(a; + 1) = yn+i{x), it follows that z{x + 1) = z{x), which completes the proof. D 



To complete the proof of Theorem |2.5[ it remains to investigate the case where 

max /3 — min /3 > Vqk sin a -\- 5 (60) 



for some S > 0. The comparison principle (Proposition 2.3) implies that stationary solutions act 
as barrier and prevents the motion of the drop down the inclined plane. However, condition ( 60 ) 
does not automatically imply the existence of a stationary solution. Nevertheless, we can show 
that stationary solutions do not exist provided the period of /3 is small enough. 



Proof of Theorem 2.5 (ii). We denote by L the period of , 



We first fix Iq such that H{£q) = max/3. We note that we can always find such an £q. 



Indeed, H{£) — > +oo as £ — >■ 0, and using the fact that G{£c) = together with (30), we get 
H{£c) = Vb'tsina. Condition ([60| thus implies 



H{£c) ~ VoKsina < inax/3 — min/3 < max/3. 

and the intermediate value theorem implies the existence of £q < £c such that H{£q) = max/?. 

We now fix a such that ^(a) = min/3, and we find h S [a+£Q,a+£Q+L\ such that /3(&) = max/3 
(we can always find such a h thanks to the periodicity of /3). 

Consider now the solution of ([7]) with support (a, h). Since £o < b — a < £q + L, we have 

i|u,|2(6) = H{b-a) < H{£o) = max;9 
and so (thanks to the choice of b) 

l\u,\\b)-m<o. (61) 



We also have (using Lemma 3.6 and the fact that £q < £c) 
1 



2 WWa) = G{b -a)> Gi£o + L) > G(4) - CL, (62) 



as long as i < \{£c — d-o)- Finally, (30) implies 



G{£t)) = H{£q) — VoKSvaa = max/3 — Vo^^sina. 



In particular ( 62 ) yields 

-T.Wx^ia) > P{a) =niin^ 
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as soon as 

max P — min P >Vqk sin a + CL 



which holds for all L < 5^ in view of ( 60 1 . Therefore we deduce 

1 



P[a)-^-\u^Wa)<G (63) 



The comparison principle together with (61) and (63) implies that if a(0) < a and 6(0) < 5, 
then the corresponding solution will satisfy 

a{t) < a and b{t) < b for all t > 0, 

hence the result. □ 

7 Homogenization of the velocity law 

We now prove our last result, Theorem |2.7[ Recall that the function r{q) is defined for all q 7^ 



as in ( 20 ) , as the effective speed of the global (periodic) solution of the ODE 

x'{t) = q-P{xit)). (64) 

(see Lemma |2.6[ whose proof is provided in appendix [A|). 

Now, let u^(x, t) solve the inhomogeneous problem with period e and with given initial data 
Mo- We also denote by a^{t) and 6^(t) the left and right cndpoints of the support of u^. 

Convergence. First, we obtain a priori estimate on m^ in the same way that we did in the 
construction of a solution in Section 14.21 

We recall that a'^{t), b^{t) and i'^{t) = b'^{t) — a'^(t) are Lipschitz continuous (bounded speed 
of the endpoints) uniformly with respect to e and there exists £ independent of e such that 

|<r(i)<4 yt>0. (65) 

In particular there exists a subsequence £ — > and some Lipschitz continuous functions a{t) 
and b{t) such that 

a^(i) — > a{t), b'^{t) -^ b{t) uniformly with respect to i € [a, b] 

for all a < 6 e M. 

Next, we note that for any t > 0, u^{x,t) solves 

—Uxxi'it) — M^^i'^)) ^ K cos au + (x — b'^{t))Ksma, in {a'^ , ¥) 

and so the function v'^{x,t) — u^{a^{t) + x£^{t),t) solves 

-Vxx{-,t) =£^tf[Xie''{t))-KCOsav + e^{t){x-l)Ksina], in (0,1) 

with 

v%0,t) ^v%l,t) =0. 
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In particular, it is readily seen that x i— ?> v^{x, t) is bounded in C^(0, 1) uniformly with respect 
to t, and that v^ and w^ are Lipschitz continuous with respect to t, uniformly with respect to e. 
Finally, v"^ converges locally uniformly (with respect to x and t) to a function v solution of 

-Vxxi-,t) ^ e{tf [X{£{t)) - KCOsav + £{t){x~ l)Ksina] , in (0,1) with v{0,t) = ^(l,^) = 0. 



Writing 



we deduce: 






Lemma 7.1. The function x i— )• u^{x, t) is bounded in C'^{a^ , b^) uniformly w.r.t. t, and u^ and 
u% are Lipschitz continuous with respect to t, uniformly with respect to e. Finally, w^ converges 
locally uniformly (with respect to x and t) to a function v solution of 

—Uxx{', t) — X{£{t)) — Kcosau + XKsina, in (a(t), b{t)) with u{a{t), t) = u{b{t), i) = 

and 

/ u{-,t)dx — Vq. 

We can now prove Theorem |2.7| 



Proof of Theorem 2.7 We first consider the equation for fc^(i), which reads 

b''{t)^H{£^t))-l3{b^t)/e). 

Since i i— ^ ^"^ is a Lipschitz continuous function, we can also write 

b''{t)=q%t)-f3{b%t)/e) 

with q'^(t) uniformly (w.r.t. e) Lipschitz function converging (uniformly) to q(t) = H{i(t)). 
For a given to > 0, the function 

x^{t) = -b'ito + et) 

solves 

x"{t) = q'ito+et)-l3{x^t)). 

Since q'^{to + et) is Lipschitz and q^{to) converges to q{to), we have that for all (5 > 0, there exists 
£0 such that 

\q'{to + et)~q{to)\<S + Ket Ve < eo, Vt. 

and so 

\q'ito+et)^q{to)\<2S Ve<£o, V|i| < -^. (66) 

Ke 

It follows that 

x"{t)<q{to) + 26-f3{x'it)) V£<eo, ^1^1 < ]^- 

and 

x"it)>qito)-26-P{x'{t)) Ve<£o, V|t| < -^. 

Ke 
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Denoting by x (respectively x) the solution of (64 1 with q = q{to) + 25 (respectively q 
q{to) ~ 26) satisfying ^(0) = a;(0) — a;^(0), we deduce 



lit) < x^t) < x{t) Ve < £o, 0<t< 



Ke' 



x{t) < x^t) < x{t) Ve<eo, 



Ke 



<i< 0. 



(67) 
(68) 



Equation (67) implies 



r{q{tn)-25)t-l<x%t)-x''{Q)<r{q{tQ) + 25)t + l Vf < £o, 0<i< 



Ke 



and so 



K' 



riq{to)-2S)t~e<b''{to + t)^b''{to)<r{q{to)+26)t + e Ve < Eq, 0<t< 
Passing to the limit e — > 0, we deduce 

r{qito) - 2S) < '('"+'\-b('o) < ^(^(^^^ _^ 2^) < t < A 

t A 

and taking the limit t — > 0^ and then (5 — > 0, we obtain (using the continuity of the function 
q'-^ r{q))- 

riqito)) < hminf ^^ + ^^ ' ^^^ < limsup ^^ + ^^ ' ^^^"^ < r(,(^o)). 



t-i-04 



t-fOj 



Furthermore, using (68 1 instead of (67), we can show that a similar inequality holds for the limit 



i — )■ 0+. We thus deduce that b is differentiable at to and satisfies 

b'ito) = r{q{h)) 
where q{to) = H{i{to)) = i|M,(&(io), to)l'- 

The equation for a^(i) is handled in a similar fashion. Recall that 

^a-'it)<G{£%t))-f3ia%t)/e), 

with equality if \u%{a^t))\ > 0. 

If lime_j.o |w^(a^(io))| > 0, then we have |u|(a^(i))| > for i in a neighborhood of ig and for 
small e, and the argument presented above applies. If lim£_j.o |M|(a^(io))| = 0, then the above 
argument applies, but we have q(to) = and so we only get 



-a''{ta + t) + a''{tQ) <r{2S)t + e Ve < eo, 0<t< 



K' 



Passing to the limit e — )■ 0, we deduce 

ajto +t) - a{to) 
t 



< r{25) < t < 



K' 
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and taking the limit ^ — > 0^ and then 5 — > 0, we obtain (using the continuity of the function 

qt-^ r{q)): 

a{to + t) - a{to) 

limsup < r(0). 

t-»-o+ t 

and a similar limit for t — > 0+ . 



This prove that u'^ converges (up to a subsequence) to a solution of (14). (15), (21 1 



To complete the proof of Theorem |2.7[ it remains to prove the uniqueness of the solution 
of the homogenized problem, which also implies the convergence of the whole sequence u^ to u. 



For that, we will show that a comparison principle similar to that of Proposition 2.3 holds (note 



that we cannot deduce this comparison principle by passing to the limit e — > in Proposition |2.3| 
since two solutions of the homogenized problem may be obtained by passing to the limit along 
different subsequences of s) . 



Using the notations of ( 29 ) we first rewrite ( [2T| ) as follows 

J min{a'(i) + r(G(<'(i))), 4 - Ut)} -= 
1 b'{t)^r{H{i{t))). 



(69) 



We now consider two solutions of the homogenized problem ui and U2 with initial data with 
support (ai(0),6i(0)) and (a2(0), 62(0)) satisfying 



ai(0)<a2(0), and 61 (0) < 62(0). 



Wc then define 



f{t) = max{ai(t) - a2{t),bi{t) - b2{t)}. 
Since the functions ai{t) and bi{t) are in M^;J^(0,oo), f{t) is in W^;J^(0,oo) and we have 



fit) 



a'lit) - a'^it) a £2(1) > £i{t) 
b',{t) - b'^it) a £2(1) < £i{t) 



When £2{t) > £i{t), we have in particular £i{t) < £c, and so (69) gives 

a[it) = -r(G(4(t))), 4(0 > -r(G(4(i))). 
We deduce 



fit) < 



r{Gi£2{t))) - r(G(£i(t))) if 4(i) > £i{t) 



r{H{£^{t)))~r{H{£2{t))) if4(0<^iW 
and the monotonicity of G, H and r implies 

fit) < 0. 

We deduce that if /(O) < 0, then f{t) < a.e. i > which implies the comparison result, and 
the uniqueness of the solution. D 
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A The function r{q): Proof of Lemma 2.6 



In this section, we prove Lemma |2.6[ We recall that /3 is a periodic function (with period 1 for 
instance) and for a given q > 0, we consider the following equation: 



x'it) = q-l3ixit)). 



(70) 



Since /3 is a Lipschitz function, (70 1 has a unique solution for any initial data x{0) = Xq and two 
solutions can never cross. 



Our first remark is that if g S [min/3,max/?], then any solution of (70 1 will be trapped in the 



sense that x{t) S [^^(O) — l,a;(0) + 1] for all t. Indeed, periodicity of /3 implies that there exists 
xi e [x{0) - 1, a::(0)) and X2 G (a;(0), a;(0) + 1] such that q - l3{xi) > and q - l3{x2) < 0. It is 
then easy to show that x{t) S [a;i,a:2] for all t. 



In that case, the effective speed of any solutions of (70) is zero. We thus have 

r{q) = 0, q e [min/3, max/3]. 



For q > max/3 and q < min/3, on the other hand, any solution of (70) will be strictly 
monotone, with \mit^±oo x{t) = ±oo. In particular, we can then show that all solutions are 
equal, up to a translation in time. Assuming that a;(0) = (without loss of generality), there 



exist a unique tj, such that x{tc) — 1. Uniqueness for (70) implies that 1 1-> x(t) is then periodic 
with period tc and that its effective speed is given by 

1 



r{q) = - 



We then have 



Proposition A.l. The function r : q i^ r{q) defined above is a non- decreasing function. Fur- 
thermore, 

(i) q I— > r{q) is locally Lipschitz in (max/3, cx)) and [0,min/3) 

(ii) If qo — max/3, then we have 

r{q) < C{q - q^f'^ for q > q^ 

if X ^^ j5 is C^ , and 

r{q)<C{-Hq-qo))-^ forq>qo 

if X i-^ 13 is only Lipschitz. In particular r is a continuous function in [0, oo) 

Proof of Proposition \A.1\ To prove the first part, we consider q2 > qi > max/3 and we set 
rj = q2 — qi. Let x{t) and y{t) be the solution of (70) with q = qi and q — q2 respectively, and 
a;(0) = y(0) = 0. Since x'{t) > qi — max/3 = J > 0, we can define 



solution of 



h'is) = 



h{s) — y o X (s) — s, 
gi - (3{x + h{x)) + rj 



< 



91 - l^ix) 
(3{x) - f3{x + h{x)) + Tj 
91 - /3(x) 

--{h{x)+Tj) 



-1 
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Since h{0) = y{x ^(0)) — = 0, the Gronwall's lemma implies 

h{l) < Ct]. (71) 



If we denote tc such that x{tc) — 1 (so that r{qi) — l/tc), then (71 1 implies 

y{tc) <1 + Crj. 
Finally, since y'{t) > S, this implies t2, such that 2/(^2) = 1 satisfies t2 > tc — Crj and so 

^■(92) < r-(gi) + Cr] 

for some C depending on qi and r{qi). 

The constant C degenerates when r{q) goes to zero, so we need different argument to prove 
(ii): Let tj = q — qg. Assume (with loss of generality) that max/? = /3(0) (and so /?'(0) =0). If /3 
is C^, then 

l3{x)>P{0)~Cx^=qo-Cx^, for aha; 

and so 

q- P{x) < rj + Cx^. 



Let now x{t) be the corresponding solution of ( 70 ). Up to a translation in time, we can always 
assume that x{Q) = 0. We then have 

x'{t) <Tj + Cx{tf 

We deduce 

.(.)<^^ iort<iCV-vr. 

In particular, 1^ defined by x{tc) — 1 satisfies 

hence the first result. 

When P is only Lipschitz, a similar argument yields 

x'{t) < ri + Cx{t) 

which gives the second result. □ 
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